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Terahertz and far-infrared electric and magnetic responses of hexagonal piezomagnetic YMnC>3 
single crystals are investigated. Antiferromagnetic resonance is observed in the spectra of magnetic 
permeability [i a [H (to) oriented within the hexagonal plane] below the Neel temperature TV. This 
excitation softens from 41 to 32 cm -1 on heating and finally disappears above Tjy. An additional 
weak and heavily-damped excitation is seen in the spectra of complex dielectric permittivity e c 
within the same frequency range. This excitation contributes to the dielectric spectra in both 
antiferromagnetic and paramagnetic phases. Its oscillator strength significantly increases on heating 
towards room temperature thus providing evidence of piezomagnetic or higher-order couplings to 
polar phonons. Other heavily-damped dielectric excitations are detected near 100 cm" 1 in the 
paramagnetic phase in both e c and e a spectra and they exhibit similar temperature behavior. These 
excitations appearing in the frequency range of magnon branches well below polar phonons could 
remind electromagnons; however, their temperature dependence is quite different. We have used 
density functional theory for calculating phonon dispersion branches in the whole Brillouin zone. A 
detailed analysis of these results and of previously published magnon dispersion branches brought 
us to the conclusion that the observed absorption bands stem from phonon-phonon and phonon- 
paramagnon differential absorption processes. The latter is enabled by a strong short-range in-plane 
spin correlations in the paramagnetic phase. 



I. INTRODUCTION 

Spin waves (magnons) in magnetically ordered mate- 
rials can be excited by the magnetic component H(u;) of 
the electromagnetic radiation, giving rise to a resonant 
dispersion of magnetic permeability in the microwave or 
terahertz (THz) frequency region. Recently, new cou- 
pled spin-lattice excitations named electromagnons have 
been discovered in multiferroics, where the magnetic or- 
der coexists with the ferroelectric one*^ Electromagnons 
are excited by the electric component E(w) of the electro- 
magnetic radiation, therefore they can be detected in the 
THz dielectric permittivity spectra. Though they were 
theoretically predicted in 1970^ the first experimental 
confirmation appeared as late as in 2006 X These exci- 
tations were mainly investigated in the rare earth (R) 
orthorhombic manganitcs i?MnC>3 and iiJM^Os (for re- 
views see e.g. Refs.0-H), and in hexaferritesii 

Multiferroics can be roughly divided into two 
groups^"— In the so-called type-I multiferroics the ferro- 
electric (FE) order takes place both above and below the 
magnetic ordering temperature and the spontaneous po- 
larization is large. However, the coupling between mag- 
netic and electric order parameters is weak. 

A general feature of type-II multifcrroic materials is 
that the ferroelectric phase is induced by magnetic order- 
ing characterized by a particular type of incommensurate 
spiral magnetic structure. In this case the magnetically- 
induced polarization is by several orders of magnitude 
smaller than in type-I multiferroics. However the cou- 
pling between electric and magnetic subsystems is large 



and giant magnctoelectric effects are observed. The 
magnon dispersion branch in the incommensurate phase 
exhibits a minimum at the wave vector q m corresponding 
to the modulation vector of the ordered spins. In con- 
trast to the magnetic resonance (magnon at q?»0) char- 
acterized by sharp spectral features, the electromagnons 
manifest themselves as very broad spectral bands be- 
cause their activation in the dielectric spectra is closely 
related to the high density of states close to the extrema 
of the magnon dispersion branches. Since the probing 
THz radiation has a long wavelength (i.e. the wave vec- 
tor q~0), the electromagnons cannot be excited by a 
resonant single-photon absorption due to the wave vec- 
tor conservation law; in this sense polar phonons should 
be involved in the interaction process. 

An experimentally observed low-frequency electro- 
magnon in type-II multiferroics was found to be related 
to the spin waves near the magnetic Brillouin zone (BZ) 
center with q=q m 4iL or to those with q=q BZE -2q m fi 2 - 
here qezE stands for the wave vector at the BZ edge. In 
both cases the low-frequency electromagnon has a similar 
frequency as the magnon with q=0 which is expected be- 
cause all these excitations are related to the same magnon 
branch. A high-frequency electromagnon corresponds to 
an excitation of the BZ-cdgc magnons (q=qBZE) which 
can induce a quasi-uniform modulation (q«0) of the lo- 
cal electric dipole moment i 11 ' 13 As for mechanisms of the 
electromagnon excitations, some researchers claim that 
the low-frequency electromagnons are activated by the in- 
verse Dzyaloshinskii-Moriya mechanism, while the high- 
frequency one by the Heiscnbcrg exchange couplingi^ 
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Other authors believe that both types of electromagnons 
can be explained by the Heisenberg exchange coupling^ 

Formerly, it was assumed that the electromagnons can 
be activated only in type-II multiferroics due to the large 
magnetoclcctric coupling. Nevertheless, electromagnons 
were recently observed also in BiFeO^r^— which is the 
most prominent type-I multifcrroic with a rather weak 
magnetoelectric coupling. In this context, THz dielectric 
spectra of multiferroics may shed new light on the nature 
of magnetoelectric coupling. 

Hexagonal manganitcs i?MnC>3 belong to the type- 
I multiferroics. In particular, the hexagonal YMnC>3 
is ferroelectric below «1250K3 and the antiferromag- 
netic (AFM) ordering sets only below T« w 70 
The magnetic symmetry is P&cmSi and therefore the 
linear magnetoclcctric coupling is forbidden. However, 
piezomagnetic and magnctoelastic couplings, and higher- 
order magnetoelectric couplings are allowed^ - — The 
piezomagnetic coupling is characterized by a bilinear 
interaction between the magnetic order parameter and 
strain, in contrast to the magnctoelastic coupling which 
is proportional to the product of squared order parame- 
ter and strain i 22 ' 25 By using the method of optical sec- 
ond harmonic generation^ the piezomagnetic coupling 
was observed owing to the interaction between AFM 
and FE domain walls in YMnQ 3 i 23 ' 24 Switching of the 
FE polarization triggers a reversal of the AFM order 
parameter i 23 ! 24 ' 27 Higher order magnetoelectric coupling 
in YMnC>3 has been observed in several works. Excep- 
tionally large atomic displacements at TV were observed 
in structural studies and they demonstrate unusually 
strong magnctoelastic coupling^ The large spin-polar- 
phonon coupling manifests itself by a decrease of the 
low-frequency permittivity 2 ^ near Tjv which is probably 
caused by anomalous hardening of several infrared-active 
phonons 3 ^. Similar phonon anomalies were observed 
near Tjv also in the Raman spectral Ultrasound mea- 
surements on a single crystal of the hexagonal YM11O3 
showed anomalous behavior of the elastic moduli Cn and 
Cqs due to a strong coupling of the lattice with the in- 
plane exchange interaction^ 

The AFM resonance in hexagonal YMnC>3 crystal was 
first observed and briefly (without any figures) reported 
in Ref. More detailed THz studies of YM11O3 ceram- 
ics were recently published in Ref. HH- The AFM res- 
onance lies near 43 cm -1 at 4K and its frequency soft- 
ens on heating towards Tjf, where it disappears! 33 ' 34 
Three magnon branches were discovered below Tjv using 
inelastic neutron scattering (INS)i 20 ' 35 ' 36 Two of them 
are degenerated near the BZ center and their frequen- 
cies correspond to the above mentioned AFM resonance. 
Moreover, a possible existence of magnons and short- 
range correlations between spins at Mn sites in param- 
agnetic phase were indicated by INSi2£~— The magnc- 
toelastic coupling manifests itself also by a strong mix- 
ing of magnons with acoustic phonons; this leads to a 
gap in the transverse acoustic (TA) phonon branch oc- 
curring at the frequencies and wave vectors where the 



uncoupled magnon and TA branches would intersect^. 
Recent polarized INS measurements revealed that the 
excitation detected at liquid helium temperatures near 
43 cm -1 has a mixed character of magnetic spin wave 
and lattice vibration^ i.e. its contribution to both the 
magnetic permeability and the dielectric permittivity is 
possible. 

The reported piezomagnetic, magnctoelastic, and 
higher-order magnetoelectric couplings in optical, acous- 
tic and mainly INS data stimulated our spectroscopic 
study of hexagonal single crystals of YMnC>3. In this 
paper, we present results on far-infrared (FIR) and THz 
polarized spectra in this material emphasizing interac- 
tion between magnetic, electric and phonon subsystems. 
We demonstrate that strongly undcrdamped AFM res- 
onance observed near w 40 cm -1 contributes only to 
the magnetic permeability spectra below Tpf. An ad- 
ditional broad and weak absorption band was observed 
in the same frequency range in the dielectric spectra both 
below and above Tjv. In contrast to electromagnons 
which are typically observed only below 50 K, the os- 
cillator strength of this excitation significantly increases 
on heating when room temperature is approached. This 
indicates that the feature must be related to the occu- 
pation number of magnons and/or phonons. Additional 
absorption band with similar temperature behavior was 
observed also near 100 cm -1 . We will show that both 
these excitations can be explained by differential multi- 
phonon and magnon-phonon processes. 



II. EXPERIMENTAL DETAILS 

The experiments were performed using a Fourier- 
transform infrared (FTIR) spectrometer Bruker IFS 113v 
and a custom-made THz time-domain spectrometer^ 
In both experiments, Optistat CF cryostats (Oxford In- 
struments) with polyethylene (FIR) or Mylar (THz) win- 
dows were used for measurements between 10 and 300 K. 
Helium-cooled bolometer operating at 1.6 K was used as 
a detector in the FTIR spectrometer. Principles of THz 
time-domain spectroscopy are explained in Ref. l42l The 
output of a femtosecond Ti: sapphire laser oscillator (Co- 
herent, Mira) excites an interdigitated photoconducting 
switch TeraSED (Giga-Optics) to generate linearly polar- 
ized broadband THz probing pulses. A gated detection 
scheme based on an electro-optic sampling with a 1 mm- 
thick [110] ZnTe crystal permits to measure the time pro- 
file of the electric field of the transmitted THz pulse (see 
Ref. for further details). 

Hexagonal YMn03 single crystals were grown by the 
floating zone method^ 3 . Two crystal plates with lateral 
dimensions of ~ 4.5 x 5 mm 2 and with the c-axis oriented 
either in-plane or out-of-plane along its normal, were 
cut and polished to obtain highly plane-parallel samples 
(within ± 1 /im) with thicknesses of 1100 and 348 (im 
for each orientation, respectively. These crystal plates 
were probed using the THz and FIR beam in all pos- 
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(a) E±c,H±c (b) E±c,H\\c (c) E\\c,H±c 
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FIG. 1: (Color online) Complex THz spectra of YMnOs taken at various temperatures. The polarization of the THz beam is 
indicated above the plots. The resonance feature near ~40 cm -1 corresponds to the doubly-degenerated AFM mode contributing 
to the magnetic permeability [i a spectra. 
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FIG. 2: (Color online) Temperature dependence of the com- 
plex permittivity e c and permeability fi a spectra calculated 
from data plotted in Fig. Q] The solid e c curves at 100 and 
140 K result from the oscillator fit. 
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FIG. 3: (Color online) Temperature dependences of parame- 
ters of the resonances observed in magnetic /j, a and dielectric 
£ c spectra. Closed circles: frequency of the AFM resonance; 
Open squares and triangles: eigen- frequency tJdieli and oscil- 
lator strength Aeu} dietl , respectively, of the mode observed in 
the dielectric spectra in Fig. [2] The dotted line shows the 
population increase of an energy level at 66 cm -1 following 
the Bose-Einstein statistics. 



III. RESULTS 



sible geometries: E(w)_l_c, H(cj)_Lc; E(tj)_l_c, H(w)||c 
and E (w)||c, H (u;)_Lc. It enabled us to get access to the 
complex spectra of the products e a ^ a , £aHc and e c /i a as 
shown in Fig. [TJa), (b), and (c), respectively. 



At low temperatures, the peak around 40 cm -1 seen 
in the spectra of e a fi a and e c fi a [Fig. QJa, c)] but not in 
those of Ea^c [Fig.lUb)] is definitely due to the AFM reso- 
nance as it contributes only to the magnetic permeability 
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fj, a . The AFM resonance vanishes above Tm ~ 70 K. The 
data shown in Fig. [TJb) allow us to assume that \i c = 1 
in the THz range. This is in agreement with the mag- 
netic order of YM11O3 in the AFM phase: the spins are 
ordered in adjacent layers in the hexagonal plane in such 
a way that the magnetic resonances are not expected to 
be excited with H||c. Based on this assumption, we are 
able to retrieve the complex values of the permeability 
Ha and of the permittivity e c (sec Fig. [2]). 

The spectra of [x a were fitted by a damped harmonic 
oscillator and the resulting AFM resonance frequency is 
plotted in Fig. [3j a strong softening is observed upon 
heating towards Tjy. Similar temperature dependence 
was briefly published earlier ; 33 ' with the magnon fre- 
quency higher by approximately 2cm~ 1 . 

Besides the sharp AFM resonance line in the low- 
temperature H a spectra one can observe a broad dielectric 
absorption band around 40 cm -1 in the e c spectra. This 
feature is detected even above TV, where its strength re- 
markably increases with temperature. The presence of 
such a resonance in e c is qualitatively expected from a 
simple comparison of the raw data in Figs. []Ja) and ( c )- 
The accessible spectral range of the THz measurements 
for our sample is limited to ~ 60 cm , therefore we have 
performed also FTIR transmission (up to 100 cm -1 ) and 
reflectivity (up to 650 cm -1 ) measurements for all polar- 
izations. 

An example of FTIR experimental transmittance and 
reflectivity spectra obtained at 120 K and their various 
fits are shown in Fig. @] Regular oscillations observed 
in the transmittance spectrum are due to Fabry-Perot 
interferences in the plane-parallel sample; a weak mini- 
mum near 40 cm -1 corresponds to the broad absorption 
band detected in the THz dielectric spectra (see Fig. [2j- 
According to Ref. as well as according to our FTIR 
reflectivity (see e.g. Fig. 0]), the lowest frequency polar 
phonons lie above 150 cm -1 in both polarized E[|c and 
E±c spectra. Nevertheless, our simultaneous fits of the 
THz complex permittivity and FTIR transmittance and 
reflectivity data reveal several additional modes below 
these phonon frequencies. The relevant spectra are plot- 
ted in Fig. [5j Besides the sharp magnon line at 40 cm -1 
three other broad modes at roughly 10, 40 and 100 cm -1 
were used in the fitting procedure in order to account 
for the measured shape of the E||c spectra at 10 K (see 
Fig. EK). The additional modes remain in the spectra 
up to room temperature and their strength increases on 
heating. Also in E_Lc polarized spectra, two broad modes 
observed near 10 and 90 cm -1 were used for the fits above 
50 K. 

The feature observed near 10 cm -1 in both polarized 
spectra could be related to low-frequency magnons 3 -^ (cf. 
the low-frequency magnon branches shown in Fig. [7|)- 
However, the sensitivity and accuracy of our THz spec- 
tra below 20 cm -1 is limited; therefore we cannot exclude 
that it is only an artifact. For this reason we will not 
speculate about the origin of this excitation. All other 
modes appearing below 150 cm -1 are clearly observed in 
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FIG. 4: (Color online) Example of the experimental FTIR 
transmittance and reflectivity spectra of 348 /J,m thick YMnOa 
crystal with polarization E||c obtained at 120 K. Dashed- 
dotted blue line: theoretical transmittance spectrum obtained 
from parameters of the FTIR reflectivity fit (without consid- 
ering modes observed by THz spectroscopy); solid red line: 
simultaneous fit of the FTIR transmission and THz spectra. 
Dashed green line is the result of a fit of the reflectivity using 
the parameters obtained from the fit of FIR and THz trans- 
mittance. One can see that the reflectivity spectrum is not 
sensitive enough to detect the weak broad modes near 40 and 
100 cm - . Oscillations in the experimental reflectivity spec- 
trum observed below 80 cm -1 are caused by the diffraction 
of FIR beam on a small sample. 



the THz and/or FTIR transmittance spectra while the 
FTIR reflectivity measurements are not sensitive enough 
to detect and resolve these weak and broad spectral fea- 
tures (see Fig. fj]). Their origin will be discussed in the 
next section. 

The temperature dependence of the sub-THz complex 
dielectric permittivity e a plotted in Fig. [5] for 20 cm -1 
exhibits a pronounced drop below TV. Such anomaly 
is a typical feature of large spin-phonon coupling which 
occurs only in hexagonal planes of YMnOs, where the 
spins are ordered. For that reason the anomaly is not 
observed in e c (T). The AFM phase transition is accom- 
panied by unusually large atomic displacements, which 
were detected by neutron diffraction^ for this reason the 
phonon frequencies change below Tn ■ The decrease in e' a 
and e" a is mainly caused by hardening of the E\ symme- 
try polar mode seen near 250 cm -1 in the IR reflectivity 
spectra with polarization ETci 3 ^ Fits of our IR reflectiv- 
ity spectra show that the mode near 250 cm -1 hardens 
from 246 cm" 1 (at 300 K) to 256 cm" 1 (atlOK) and 
therefore its dielectric contribution Ae is reduced from 
9.1 (300 K) to 7.6 (10 K). This decrease of Ae is mainly 
responsible for the change of the permittivity £„(T) seen 
in Fig. [6l Hardening of other modes brings a minor con- 
tribution to the decrease of £„(T) on cooling. Similar 
temperature dependence of e' a was observed also in the 
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FIG. 5: (Color online) The measured THz loss spectra of 
YMn03 (symbols) and those obtained from the fits of FTIR 
transmittance and reflectivity spectra. Below Tat=70K, the 
spectra correspond to imaginary parts of the permittivity- 
permeability product. Above TV, the spectra correspond 
to the dielectric losses. Polarizations of electric and mag- 
netic components of IR or THz beams are indicated. Dashed 
lines are the fits of room-temperature FTIR reflectivity spec- 
tra without taking into account the IR and THz transmit- 
tance spectra. The marked peaks above 150 cm" 1 are due to 
phonons; the origin of lower frequency absorption bands is 
discussed in the text. 



radio-frequency region^ providing evidence of the ab- 
sence of dielectric dispersion below 100 GHz. Gradual 
decrease of e' a and e' c on cooling from 300 to 100 K is a 
usual behavior caused by a small phonon stiffening as a 
consequence of thermal contraction. 



IV. DISCUSSION 

The question arises about the origin of the absorp- 
tion bands appearing below phonon resonances in Fig. [5] 
They are much weaker and significantly broader than 
those of polar phonons and their strength increases when 
the temperature is increased, i.e., the strength is high in 
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FIG. 6: (Color online) Temperature dependence of the (a) 
permittivity and (b) dielectric loss measured at 20 cm -1 with 
polarization E±c (red solid lines) and E||c (black dashed 
lines). 



the paramagnetic phase. Their frequencies lying in the 
range of 40-100 cm" 1 coincide with those of the magnon 
branch observed by INS at 7 K over the BZ2£ (see Fig.©. 
In the following text we discuss whether these features 
can be related to the magnon dispersion branches. 

Could a spin wave still exist in hexagonal YMn03 at 
room temperature? It is well established that Mn spins 
exhibit a strong short-range correlation in hexagonal 
YMnOa far above Tn- This was proved by an anomalous 
behavior of the thermal conductivity^ elastic moduli^ 
as well as by neutron scattering experiments Nev- 
ertheless, due to the short-range correlation of the spins 
in the hexagonal plane of YM11O3, one can expect the 
existence of only short- wavelength paramagnons, i.e. 
magnons with large wavevectors q x near the M-point of 
the BZ. A part of such a paramagnon branch is schemat- 
ically plotted in Fig. [7] Note that its frequency is lower 
than that of the magnon branch at 7 K, as the magnon 
frequency decreases by almost 10 cm -1 on heating to- 
wards T w (see Fig.©. 

Electromagnons are excitations with frequencies close 
to those of spin waves, which, due to specific couplings, 
are activated in the dielectric spectra. In pcrovskite 
manganites, the parts of magnon branches exhibiting a 
high density of states are mainly involved in these in- 
teractions (at BZ edge or close to the spin modulation 
wave vector)^ However, these electromagnons were ob- 
served only at very low temperatures (typically less than 
50 K). Their strength dramatically decreases on heat- 
ing and they usually disappear from the spectra at Tjv 
or close above Tjv^— This is in contradiction with our 
observations in YM11O3. 

We came to the conclusion that the broad absorption 
bands we observe in the dielectric spectra reflect excita- 
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FIG. 7: (Color online) Dispersion branches of phonons (the- 
oretical; black solid lines) and magnons (experimental at 
7K; red dashed lines). The red-dotted line indicates the 
presumable dispersion of the paramagnon near the M-point. 
The symbols shown at the BZ edges indicate the polariza- 
tion of the phonons at the BZ boundary: a and c stand for 
phonons polarized within the hexagonal plane and in the per- 
pendicular direction, respectively. In the T-point, the Ei- and 
4i-phonons observed experimentally^ are marked by green 
and blue points, respectively; other modes are silent. Blue 
arrows with assignment uodieii and uidiei2 indicate phonon- 
paramagnon excitations observed in the dielectric loss spectra 
of e". Green arrow marked as LJdieis indicates a broad multi- 
phonon absorption observed in the e" loss spectra (see Fig. [5]). 



tions which must be coupled to phonons. Let us discuss 
in brief which types of interaction between the magnetic 
subsystem and other degrees of freedom might be ex- 
pected on the basis of the point group crystallographic 
symmetry 6mm and the magnetic symmetry 6mm..— The 
magnetic order parameter of YMnC>3 was analyzed in 
several publications and it was shown to transform fol- 
lowing l$i (T4) irreducible representation of the 6mm 
group 
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The 6mm symmetry strictly forbids the lin- 
ear magnetoelcctric effect, i.e. bilinear terms onjHiEj 
(where Hi and Ei arc components of the magnetic and 
electric field, respectively) are not allowed in the ther- 
modynamic potential^ However, a higher order magne- 
toelectric effect (called sometimes the magnetodielectric 
effect), accounted for by the PijkHiHjEk terms in the 
thermodynamic potential, is allowed. This effect mani- 
fests itself in our measurements as a kink near Tjy in the 
temperature dependence of e' a (see Fig. IB]). 



The magnetic symmetry of YMn03 allows the piczo- 
magnetic contribution to the thermodynamic potential 
described by the terms PijkHitjjk, where ajk is a stress 
component and pijk denotes components of the piezo- 
magnetic tensor. 22 ' 24 We believe that this type of bilin- 
ear coupling must play an important role in the inter- 
action between the magnetic subsystem and the lattice. 
Usually, the piczomagnetic effect is allowed thanks to 
the relativistic part of spin-lattice and spin-spin inter- 
actions, provided the symmetry restrictions are met.— 
However, in YMnC>3 which is a noncollinear antifcrro- 
magnet the exchange (Coulomb) interactions may be 
by several orders of magnitude stronger than the rel- 
ativistic ones and, therefore, they can be the origin 
of piezomagnetism.— For example, extraordinary spin- 
phonon interactions were shown to contribute to the ther- 
mal conductivity of YMnO-3 below Tat.— Higher order 
effects such as pijkiHiHjaki are naturally also allowed in 
YMn0 3 . 

In order to provide a more quantitative explana- 
tion of the interaction between magnetic subsystem 
and phonons, we calculated the phonon spectrum from 
first principles within the spin-polarized local density 
approximation^. We used projector augmented-wave 
potentials as implemented in Vienna Ab Initio Sim- 
ulation Package (VASP)^r— . The following valence- 
electron configurations were considered: 4s 2 4p 6 5s 2 4d 1 
for Y, 3p 6 4s 2 3d 5 for Mn, and 2s 2 2p 4 for oxygen. To ac- 
count for the strong electron correlation effects on the 
rf-shells of Mn atoms, we used LDA+U approac h 56 i 57 
with an on-site Coulomb parameter U = 8.0 eV and 
Hund's exchange Jh = 0.88 eV as calculated in Ref. l58l . 
The spin-orbit interaction was not taken into account. 
We used an A-type antiferromagnetic structure, where 
spins on Mn honeycomb layers are aligned ferromag- 
netically and the layers with opposite spin-direction al- 
ternate along the c-axisi 50 ' 59 A kinetic energy cutoff of 
500 eV and a 4 x 4 x 2 T-centered fc-point mesh was 
used in the structural relaxation of the unit cell, where 
the Hellman-Feynman forces were minimized to a value 
smaller than 0.5 meV/A. Phonon calculations were per- 
formed on a 2 x 2 x 1 T- centered fc-point mesh, with a 
2x2x2 superccll within the force-constant method i 60 ' 61 
The Hellman-Feynman forces were calculated for dis- 
placements of atoms of up to 0.04 A. The dynamical 
matrix for each g-point in the BZ was constructed by a 
Fourier transformation of the force constants, calculated 
for the T-point and for the BZ boundaries. Phonon-mode 
frequencies and atomic displacement patterns for each q- 
point were obtained as eigenvalues and eigenvectors of 
the dynamical matrices. The result for directions A-r-M 
and wavenumbers up to 200 cm" 1 are presented in Fig. [7] 

As we have already pointed out, the absorption 
strength significantly increases on heating. This is typi- 
cal for difference frequency absorption. Such process in- 
cludes the annihilation of one quasi-particle (phonon or 
magnon) with frequency lu\ and the creation of another 
quasi-particle with a higher frequency uji- The dielec- 
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trie resonance then occurs at frequency LOdiei — ^2 — ^1 ■ 
This process can involve excitations from the whole BZ 
provided that the total wave vector is conserved. The 
contribution of the parts of the dispersion branch with 
the highest density of states is expected to dominate. 
The high number of available states is found namely at 
the flat parts of the bands close to the BZ boundaries, as 
it was observed, for example, in MgO£2- 

Obviously such a process is strongly temperature de- 
pendent, as it is related to the population of excitations 
with frequency uj\, which follows the Bose-Einstein statis- 
tics. At low temperatures the population of the lev- 
els which we study is close to zero and the differential 
absorption then practically vanishes. It becomes more 
probable when the energy level is thermally populated 
at higher temperatures. This is in qualitative agreement 
with our observations. 

The differential transitions at the BZ boundary are 
possible only between phonons with the same symme- 
try and if the total wave vector is conserved (i.e. the 
transition must be vertical in the wave- vector space). 
The broad absorption around 90 cm -1 seen in e a spectra 
(Fig. [SJd) can be explained by differential multiphonon 
absorption. Phonons near 60 and 150 cm -1 at the A- 
point of BZ are polarized in the hexagonal plane (marked 
as a in Fig. [7]) and their difference gives the frequency 
= 90 cm -1 , as observed. 

However, the two bands seen in e c spectra around 
Wdieii — 40 cm -1 and u>diei2 = 100 cm -1 are impossi- 
ble to explain by multiphonon absorption. The fre- 
quency of the c-polarizcd phonons at the BZ edge is 
higher than 100 cm -1 . It means that the population of 
such phonons should be much lower than that of the 
a-polarized phonon at 60 cm -1 . For this reason the 
strength of the differential multiphonon absorption in the 
e c spectra should be weaker than in e a spectra. Moreover, 
within such a hypothesis, a continuous absorption band 
would be expected in the spectra due to the large num- 
ber of c-polarizcd phonons at the M-point (see scheme in 
Fig. [Z])- This is in contradiction with the experimental 
results presented in Fig. [SJ 

We assume the existence of paramagnons near the M- 
point, and in this case a differential paramagnon-phonon 
absorption with several maxima can be obtained. More- 
over, because of the similar Bose-Einstein factor for the 
paramagnon close to 70 cm -1 and phonon near 60 cm -1 
at the A-point, the absorptions observed in both e a and 
e c should have comparable strengths. This fits well with 
the experiment. The frequency uJdieli increases on heat- 
ing (Fig. [3]) presumably due to the softening of the para- 
magnon branch with increasing temperature. The in- 
crease of the oscillator strength Asujj iell of the mode 
observed in Fig. 3 is compatible with the temperature 
increase of the Bose-Einstein factor: this is demonstrated 
by the dotted line which shows the expected population 



increase of an energy level at 66 cm 1 (i.e. the frequency 
of paramagnon at qezE)- 

V. CONCLUSIONS 

The THz and FTIR transmission spectra of hexagonal 
YMnC>3 clearly revealed two kinds of excitations of differ- 
ent nature, which exist below polar phonon frequencies. 
The sharp AFM resonance band observed near 40 cm -1 
at low temperatures broadens upon heating and disap- 
pears close to T^r. This resonance is the main contrib- 
utor to the magnetic permeability fi a . Additional broad 
excitations were observed in the frequency range 40-100 
cm -1 in the dielectric permittivity spectra in both the 
AFM and paramagnetic phases. Our theoretical expla- 
nation of the activation of these excitations in the THz 
dielectric spectra is based on a two-particle differential 
processes schematically shown in Fig. [7J The resonance 
observed in e a spectra is caused by differential phonon 
absorption in the A-point of the BZ. The two broad ab- 
sorption bands in e c spectra were described as differential 
phonon-paramagnon processes. The absorption strength 
of these excitations in the THz spectra increases on heat- 
ing due to the growing population of paramagnons and 
phonons with temperature. This is possible in the para- 
magnetic phase owing to strong short-range spin corre- 
lations within hexagonal planes of YM11O3. The pro- 
cesses we observe in YMn03, where the linear magneto- 
electric coupling is forbidden, are clearly different from 
the one responsible for the appearance of electromagnons 
in multiferroics with spin-induced ferroelectricity.-^>^ 
The multiphonon absorptions are allowed by symmetry 
in all dielectric systems, while paramagnon-phonon ab- 
sorptions can be expected only in paramagnetic systems 
with a strong short-range magnetic order (e.g. in hexag- 
onal manganites). Magnon-phonon absorption should 
be also detectable in all magnetically ordered systems 
(FM, AFM, ferrimagnets etc.) with relatively high criti- 
cal temperatures. In such conditions the magnons at the 
Brillouin zone edge may become sufficiently populated to 
allow multiparticle effects in the spectra. This may stim- 
ulate further THz and FIR studies of other magnetically 
polarizablc systems. 
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